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ABS TRACT 


The solution of a beam on an elastic nonlinear foundation by the 
Finite Element Method is presented in this thesis. Discontinuous 
(Winkler) and continuous foundations are considered. The general 
formulation is based on the Galerkin method. The solution technique 
for the linear case uses Gauss elimination and the solution technique 
for the nonlinear case is based on Brown's method (a modified Newton- 


Raphson). Some illustrative examples are presented. A brief com- 


parison of continuous and non-continuous foundations is made. 
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I. INTRODUCTION 


In treating the various problems related to structures on elastic 
foundations, it has been customary for engineers to assume that the 
pressure in the foundation is linearly proportional at every point to the 
deflection of the foundation at that point, independent of the pressure or 
deflection occurring in other parts of the foundation. This assumption, 
is mathematically by far the simplest that one can make regarding the 
nature of a supporting elastic medium. It assumes a complete lack of 
continuity in the material of the foundation as if it consisted of a series 
of independent springs which deflect when directly loaded. This theory 
has proved adequate for calculating stresses and deflections in railroad 
tracks and has found many applications to plate and shell structures 
[Ref. 1]. However no particular claim has been made that the deforma- 
tion or pressure distribution in sree earth foundations could be pre- 
dicted by this method. 

In recent years, some authors have treated foundations as a con- 
tinuous medium which, in contrast to the Winkler hypothesis, repre- 
sents the case of complete continuity in the material [Ref. 9]. It is 
evident that the two types of foundation models lead to different results, 
but how quantitatively significant is the difference has not previously 


been made clear. 


iz 





The properties of actual foundations, however, seem to lie some- 
where in the gap between these two extreme cases. The physical prop- 
erties of soils are obviously of a very complicated nature. Their de- 
formation behavior is influenced by a number of factors such as the 
physical structure, porosity, existence and movement of fluids in the 
pores, etc. In addition, such geologic features as faults, joints, seams, 
crushed zones, fissures and other tectomic effects produce behavior 
significantly different from that derived on the assumption of continuous 
mass [Ref. 14]. The question is then raised: does the approximation 
based on the Winkler hypothesis still hold over various types of soil 
foundations? If it does not, then how much error may be expected? 
How much does the shearing effect in the material of the foundation 
contribute to its deformation? Considering only the deflections of one- 
dimensional beams on foundations (for simplicity), this thesis attempts 


to answer these questions. 


13 





II. DESCRIPTION OF ELASTIC FOUNDATIONS 


This chapter presents brief descriptions of the Winkler, Modified 
Winkler and Continuous foundations considered in this thesis. Many 


other foundations have been proposed over the years by numerous 


investigators. Kerr [Ref. 31] presents a summary of a number of 


linear foundations. 
A. DEFINITION OF WINKLER, MODIFIED WINKLER AND 

CONTINOUS FOUNDATIONS 

The following definitions will be used throughout this thesis. 

1. Classical Winkler Foundation 

This type of foundation is characterized by the fact that the 

deflection at every point of the foundation is linearly proportional to the 
pressure applied at that point and is independent of pressures or deflec- 
tions acting elsewhere in the foundation. This assumption is equivalent 
to considering the foundation as a discontinuous medium composed of 
independent elastic springs. It is believed that this hypothesis was 
proposed in 1867 by ape [Refs. 27, 9, 29]. 

@. Modified Winkler Foundation : i 


This type of foundation, still a discontinuous medium, is 


eSome authors claim that Euler seems to have been the first to 
formulate the hypothesis, although it is generally attributed to Winkler. 
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however characterized by the fact that the pressure at every point is 
nonlinearly proportional to the deflection at that point. 

In this thesis, reference to a Winkler foundation means either 
the classical (linear) Winkler foundation or the modified (nonlinear) 
Winkler foundation. 

3. Continuous Foundation 

In contrast to the Winkler foundation, the pressure and deflec- 
tion at a point in a continuous foundation is affected by the behavior of 
the entire foundation [Refs. 1, 9, 27]. 

4. Mixed Mode Foundation 

The real foundations of natural soils are more complex than 
either the Winkler or continuous foundations. They are neither com- 
pletely continuous nor discontinuous [Refs. 14, 6]. The pressure and 
deflection at any point are nonlinearly dependent. This type of founda- 
tion is therefore considered as a combination of the modified Winkler 


foundation and the continuous foundation. 


Eee MATHEMATICAL MODELS 

In this chapter the governing equations of the foundations considered 
in this thesis are developed. - 

Consider an elastic foundation which is a compressible single layer 
of thickness H placed ona rigid base. It will be assumed that the thick- 
ness of the elastic foundation, its support conditions, the elastic con- 


straints and all other properties as well as the external load do not 
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vary in the z-direction. A narrow plate of thickness w is cut from the 
elastic foundation by two planes parallel to the x y plane. A beam of 


the same thickness w rests on the surface of this elastic foundation. 


lateral loading 
q(x) 


{le 


Beam 


elastic foundation 





Rigid base 
Figure 1. Geometry of an elastically supported beam 
Let an external load q(x), pounds per unit length, act on the beam 
(Figure l). If r(x) is the reaction of the elastic foundation against the 


beam, the differential equation of bending for the beam is then [Ref. 8]: 


(EIV") = q(x) - x (x) (1) 
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where EI (x) and V (x) are the flexural rigidity and the deflection of 
the beam, respectively. 

Equation (1) is true for any type of foundation. It contains two 
unknown functions, V (x) and r (x). In order to determine them, an 
additional relationship is needed. This relationship is associated with 
the response of the elastic foundation and depends on the type of 
foundation being considered. 

1. Winkler Foundations 

a. Governing Equation 
The Winkler foundation is characterized by a foundation 
modulus, generally denoted by the letter k* Let V, inches and r, pounds 
per inch, be the deflection and the reaction of the foundation, respectively. 
The Winkler foundation is characterized by the equation 
r (x) = KVP(x) (2) 
For the classical Winkler foundations, p= 1 and kis a 
constant and is expressed in pounds per cubic inch. This assumption 
has proved adequate for calculating stresses and deflection in railroad 
tracks [Ref. 1]. 
b. Values of Foundation Modulii 
For sand, k*may vary from 25 to 100. lb/cu.in. For 
ordinary soils on which railroad tracks have been built, k* may vary 
from 110. to 130. lb/cu.in. For gravel, the value of is even more 


uncertain and it can vary from 200. to 1200. lb/cu.in. [Refs. 1, 27]. 
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2. Modified Winkler Foundations 
For actual foundations of natural soils, k*is the coefficient 
associated with the nonlinear term [Ref. 6]. Therefore, to generalize 
the Winkler hypothesis, let p be any positive real number. Such 
foundations are called modified Winkler foundations. 
3. Continuous Foundations 
a. Governing Equation 
Let s be the shear force in the foundation, pounds per 
inch, and be the foundation modulus, pounds per cubic inch. The 
shear force s is a measure of the friction force between soil particles. 
The governing equation of continuous foundations is [Ref. 9]: 
r (x) = KV (x) - (s V (x))' (3) 
Details of the continuous foundation analysis are given in 
Appendix A. 
b. Values of s and k 
The values of s and k depend on the contact area of beam 
and foundation per unit length of beam. Fora loaded area of 6.5 inches 
of width, equal to that of a 12WF27 standard beam, and 1. inch of length, 
typical values of s and k for some soils are: 


Average sea floor sediment: 


oS 18. Ge a, 2 ohoy ies 
A typical beach soil: 
s = 1000. Ib [a=t5008 Iby inc 


where k = wk” 
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A typical inland soil: 

s = 6000. lb k = 3000. iene 

Additional values of soil characteristics are given in 
Reference 7. 

4. Mixed Mode Foundations 
a. Governing Equation 

The governing equation of this type of foundations can be 
derived from the definition of the mixed mode foundations given in 
part A. The mixed mode foundation is a combination of a modified 
Winkler foundation and a continuous foundation: 

' ! 

r (x) =k VY (x) - (s V (x) (4) 

It is seen that any of the previous foundations are special 
cases of the mixed mode foundation. For the Winkler foundation s = 0 


and p = 1, for the nonlinear foundation s = 0 and for the continuous 


foundation P = l. 
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lil, GENERAL FORMULATION OF THE PROBLEM 


A. GOVERNING EQUATIONS 
Consider a beam on an elastic foundation. The differential 
equation of bending for a beam on an elastic foundation is given by 
[Ref. 8]: 
(EIV"(x)) = q(x) - r(x) (1) 
where q(x) is the external loading, V(x) is the beam deflection and r(x) 
is the reaction of the foundation. The mixed mode foundation, by 
definition, can be considered as the most general type of foundation. 
Its governing equation is given Eqn. 1 with 
-~,vP ey) 
r(x) = kV% (x) - (sV (x)) (4) 
and p> 0. 
Substitution of (4) into (1) gives: 


(EIV (x)) = q(x) - kV? (x) + (sV (x)) (5) 


" ! 


(EIV") - (sV') +kvP =q (6) 


To solve this equation for any positive real p , consider a least 


square best fit polynomial such that 


ees Sac bo NAGS Tec b, V(x) te een (7) 


where the b. coefficients (i = 1,2,3,...) are constants. 
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For simplicity in the development of an actual computer program, 
equation (7) is taken to be a second degree polynomial. This restriction 
is not inherent in the Galerkin procedure but rather in the Finite Element 
formulation of Galerkin adopted in this thesis. Substitution of (7) into 


(6) leads to the following equation: 


(EIV")" - (sV')' +k (by + bZV +b,V%) =q (8) 

Or 
Tt 1 

(EIV") - (sV') +kboV 4+ kb, V° = al < ay (9) 
Letting 

a, = kb. (10) 

a, = kb, (11) 

f = q(x) - kby (12) 
equation (9) may be rewritten: 

2 
(EIV") - (sV')' +a,V 4 2p a (13) 


1. Beam ona Winkler Foundation 
Equation (13) can be applied to the case of a beam ona classical 


Winkler foundation, by letting bj] = b, = 0, bz = 1 andS), =0. Hence 


3 
the governing equation of bending for a beam ona classical Winkler 
foundation is: 
(EIV")) +a,V = q(x) (14) 
where a, is defined by Eqn. 10. 
2. Beam ona Modified Winkler Foundation 


When s = 0, Equation (13) becomes the governing equation 


of bending for a beam ona modified Winkler foundation: 


eal 





(EIV") +a,V + are = f (15) 
where terms are defined by Eqns. 10, 11, ld. 
3. Beam ona Continuous Foundation 


Equation (13) can be applied to the case of bending of a beam 


=e ieee. = Ole: 


ona continuous foundation, by letting b ] 3 


Z 
(EIV") - (sV')' +a7V = q(x) (16) 
4. Beam ona Mixed Mode Foundation 
Equation (13) is the governing equation of bending for a beam 
On a mixed mode foundation: 


(EIV") - (sV') +a,V + qo V“ yy; (17) 


B. BOUNDARY CONDITIONS 
Reference 8 thoroughly develops the boundary conditions for a beam. 

Only the boundary conditions on displacement and slope (the so called 

essential or principal boundary conditions) must be imposed ina 


variational formulation of the problem. 


e ASSUMPTIONS AND RESTRICTIONS 

Some restrictions have already been stated in the introduction 
section of this thesis: other restrictions will now be discussed. 

1. The flexural rigidity EI must be "slowly varying" in accord- 
ance with the development of the basic equations of classical beam 


theory. 


Ze 





2. The shear coefficient s must also be "slowly varying" in 
accordance with the development in Appendix A. 
: 3. The foundation modulus k(x} and the load variable q(x) need not 


be slowly varying. 


D. SOLUTION TECHNIQUES 
The Finite Element Method via the Galerkin approach will be used 
here. Its formulation will be presented in the next section. There are 
two types of problems to be considered. 
1. Linear problems 
Consider the linear system 
fama) visi Hew = 4 , (18) 
Suppose the solution for q = q) is V, and the solution for q = 
qd2 is V5, then forq=q,+ @ qd, we have the solution: 
Vee Vga cv we), 
The application of the Finite Element Method to any linear 
differential equation leads to a linear system of algebraic equations 
which will be solved by Gaussian elimination. 
2. Nonlinear problems —— 
In contrast to the linear case, any nonlinear differential 
equation, like (13) or (15), leads to a nonlinear system of algebraic 


equations which, in this thesis, will be solved by Brown's iteration 


as 





method. Reference 13 thoroughly develops the Brown's iteration method 
which is a modified Newton-Raphson method. 

The results of the Finite Element Method solutions presented 
here will be checked with the classical solutions if available, or with 
other methods, suchas the Finite Difference Method. Programming 


details will be presented in Chapter V. 
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IV. BINITE ELEMENT PORrMuUbA TION 


mo APPROXIMATION BASED ON THE GALERKIN PROCEDURE 
In order to obtain a numerical] solution for displacements of beams 
on elastic foundations, let V be approximated by the m degrees of 
freedom expression: 
Vx » V. G; (x) (21) 
= Ih 
where G;(x) are global (or system) shape functions and the unknowns 
WS coefficients are generalized coordinates [Ref. 14], i.e., system 
degrees of freedom. Since the Galerkin method deals directly with the 
differential equation [Refs. 21, 15, 22], (21) will be inserted in (13) 
and the equation residual is formed as follows: 
ae ats 
py |(ENVG mene) > = Ce VG) Gen) enc 
i= 1 
m 
+ Cay ViGAx), Geix) > + Dy <a 
mf, G(x) > =0 


IG, 
where the notation < a, b> means fatxarb(xJax . Then 


m t 
2 LV, < (EIGY 0), GG) > - Ve < (8 Gi GO), G(x) > 


mM 
ta,V, < Glx), Glx)> +4 2 a5 VV < GiesG,(x), Suey 


eral, G,.(x) » =0 


os 





Integration of the first two terms by parts yields: 


se 1" " 
rb [V. < EIG (x), G(x) > + Vi sGilx), G(x) > 
7 m 
$asV, < G(x), Gix)> + 2, agViV, < Gilx)G;(x), Gtx) >] 
- < £, G(x)> =0 (22) 


There is a constant of integration due to boundary conditions left 
by the process of integration by parts in Eqn. 22. However this constant 


can be omitted in the Galerkin procedure [Ref. 23]. 


Let 
vy if 
Sa = ¢ EIG. (x), G(x) > 
I ! 
Ci, - <sG.(x), G(x) > 
(23) 
MC < G(x), G(x) > 
Then, Eqn. 22 may be rewritten as follows: 
m 
i=] 
ra m | 
t+ 2 Bo 8, Nix ViVj- % = 9 — (24) 
i=] j=l 
foe = 1,2,3, ..., Mm 


This is a system of nonlinear algebraic equations, where the 
unknown displacements V; must also satisfy the imposed boundary 


conditions. 


26 





Be ELEMENT SHAPE FUNCTIONS 
Consider a beam element. There are two nodal coordinates, the 
deflection and slope at each end, and therefore each element must have 
a total of four degrees of freedom. On the element level these coor- 
dinates are numbered as follows: 
Coordinate 1 = V{(o) Coordinate 3 = V(t) 
Coordinate 2 = V (0) Coordinate 4 = V (4) 
The compatible element shape functions must be cubics and it is 
evident that four independent shape functions are required [Refs. 14, 16]. 
Let G. (x), the element shape functions, be defined as follows, where the 


superscript "e'' refers to the beam element: 








t 
e Z 3 
V(o) G, =l.- ane + 236 
, i+ +3 
42 i 
: x 
7 wn Gy = er 
32 3 
‘en = x? = i 
t : Rarer, 
—— } fT 


vi) 


Figure 2, Element shape functions 


e 
Note that the function G. (x) represents the variation of V along the 


element in such a way that G{(o) = G$(o) = G&(1) = G(r) = 1 and that 
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: € e e! 
other evaluations of G (0), G; (1) a G 
i 


(Ge oy (1) vanish; cf. the sketches 
i 


in figure 2. where x is the local coordinate from left to right and 1 


denotes the length of the beam element. 


C. FORMATION OF ELEMENT MATRICES 


Once the element shape functions are defined, the components of 


c 


the element matrices Kae one M- ; On and Nise can be computed 


k ij 
through ordinary integration of polynomials. The following results, 


associated with the linear terms, have been given in many text books 


vets. 21, 14, 25]. 


i eS We ode 

; fae 

K° = <EIG® (x), G® (x)>= EIS fd ie oe 
ik 1 “Sic } 

Symmetric Ze 208 

So Vc 


cae 


where EI is a constant (the average EI) over each element. 


ee Oot i is 0.1 
i i 
e e! e! e 7 : 
Ci =<8 G,(),G, (x)= s . 1333331 -0.1 -0.033333% 
symmetric eZ 2621 
r 
0.133333 } 


e . 
where S_ is a constant (the average) over each element. 
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Ms Z 
ST14297 evo zoe lt . 128571) 2 0e0e52r 


e e e 3 Z 
Mi =  < Gi (x), G(x)> = .009524t .030952F  -.0071431 
° 2 
.371429t -. 0523811 
Symmetric 
3 
. 0095244 
54 
Q® = Ga = f 0833331 
= <br OD = 7 
5h 
2 
-. 0833333 


e 
The components of eile associated with the nonlinear response are 


not inthe literature. They are as follows: 


e € e e 
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e 
N44 


e 
Notice that, according to Eqn. 23, the N.. 
/ 


: 003175t 


: 00119114" 


3 
mOOSTT St 


4 
-. 0007941 


2 
.0170641 


_. 0035722" 


& 
. 000794} 


. 3071431 


ae 0384921° 


: 006349t- 


4 
= -. 0011901 


‘ie coefficients are equal 


for any permutation of i, j, Kk (i.e., N112 = Ni2] = N511> etc. ) 


D. 


FORMATION OF SYSTEM MATRICES 


The system matrices are obtained from the element matrices as 


follows: 


A correspondence table between local and global coordinates 
is formed 


The element coefficients of any matrix, say the element 
e 


stiffness coefficients Kip 


are substituted directly into the 
system stiffness matrix Ks according to the correspondence 


of element coordinates i, j to system coordinates I, J.[Ref.14] 
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Ve. COMPUTER PROGRAMS 


Equation (13) is originally of the form: 

(EIV") = (sV'). + kv? =q (29) 

The following block diagram shows how the investigation proceeds 
(Figure 3). Equation 13 or 29 will be solved by the Finite Element 
Method. First consider the case of a classical Winkler foundation, 
for which the theoretical analysis is readily available. The results 
from the Finite Element Method will then be compared to the theoretical 
results. Similarly, the modified Winkler foundation will be solved, also 
by the Finite Element Method. Since there is no theoretical analysis 
for the case of these nonlinear problems, some of the results from the 
Finite Element Method will be compared to those from the Finite Dif- 
ference Method. Finally, in dealing with continuous foundations, cer- 
tain types of soils will be considered. Results obtained by varying 


flexural rigidity of beams, foundation modulus and applied force will be 


presented in the next chapter. 


A. GENERAL PROGRAM ORGANIZATION ‘ 

The program may be divided into two separate categories, depend- 
ing on linearity or nonlinearity. The reason for this distinction is that 
for nonlinear problems, the Brown's iteration technique is applied and 


double precision is required. In contrast, a simple Gaussian elimina- 


tion method may be applied to solve the linear systems. 
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Figure 3. Sequence of the problem investigation 
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Figure 5. Matrix formation (STIFF1 is used for classical 


Winkler foundation only, STIFF2 is used for all other cases. ) 
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Figure 6. Actual problem solution: 
(a) for linear problem, 


(b) for nonlinear problem. 
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Figure 7. Finite Difference Method solution 
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Each type of program performs four major functions: 
1. Storing all information needed for subsequent problems 
2. Supplying necessary data to each particular problem 
3.a. Forming the element matrices 
b. Forming the system matrices 
4. Solving the algebraic problem. 
This organization is shown in Figures 4 through 6. 
The Finite Difference Method is applied to some problems to thesk 


the results. Its general flow chart is shown in Figure 7. 


B. SUBROUTINES 
The MAIN program is a control segment that serves to call other 
subroutines in proper order as required. No calculations are performed 
in the MAIN program (See Appendix D). 
1. Subroutine STORE 
This subroutine reads all data available for the problem under 
investigation, including the maximum number of elements, the number 
of nodes, the number of boundary conditions, the number of degrees of 
freedom, the correspondence table of global and local nodes, load 
conditions, beam properties and foundation characteristics. 
2. Subroutine TEST 
This subroutine is used particularly for testing the convergence 
of the solution. In doing so, it selects special information and supplies 
it to the MAIN program which, keeping all information constant, doubles 


the number of elements from one run to another. 
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3. Subroutine .G Uke iL 
For the case of a nonlinear elastic foundation for which the 
reaction r(x) = Kye , where Pp is not an integer, a direct application 
of the finite element method via the Galerkin procedure is not efficient. 
Using a least square best fit, this subroutine replaces kvP by a 
second order polynomial. 
4, Subroutine INCHK 
Depending on various conditions in changing beam flexural 
rigidity and foundation modulus, this subroutine selects information 
from subroutine STORE and supplies it to the MAIN program. This 
subroutine also prints the appropriate echo check. 
5. Subroutines STIFF1 and STIFF2 
Subroutine STIFF1 is used for the linear problem and uses 
Single precision; subroutine STIFF2 is used for the nonlinear problem 
and uses double precision. Both of these subroutines form the element 
stiffness matrices and assembles the system stiffness matrix. 
6. Subroutine LOAD 
This subroutine forms the element load vector and assembles 
the system load vector. It can be used for both linear or nonlinear 
problems, provided it is in double precision for am later case. 
7. Subroutines BOUND1 and BOUND2 
These subroutines apply the boundary conditions to the system 


force vector and the system stiffness matrix produced by subroutines 


38 





LOAD, STIFF1 or STIFF2, and transform them into an appropriate 
form ready to be solved. Subroutine EOUND1 is used for linear problems 
and subroutine BOUND2Z for nonlinear problems. 
8. Subroutine SOLVE 

This subroutine employs the Gaussian elimination method to 

solve the linear system of algebraic equations for the linear problem. 
9: Subroutine RESULT 

This subroutine prints any output data from either linear or 
nonlinear problems. For nonlinear problems, double precision must 
be specified. 

10. Subroutine DIFSOL 

This subroutine solves the nonlinear problem by the Finite 
Difference Method of which the development will be presented in 
Chapter V, part D. This solution is optional and was employed ina 
few cases in order to check the results of the nonlinear problem given 
by the Finite Element Method. 

11. Subroutine FSOIL 

Similar to subroutine CURFIT, using a least square fit, this 
subroutine replaces the total settlement curve of a natural soil given 
by experimental data [Ref. 30], by a polynomial of second degree. The 
total settlement curve of the sea floor sediment is supposed to have the 
same shape as that of the inland soil but its magnitude is reduced by a 


certain scale. 
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12. Subroutine VGUESS 
This subroutine is designed to give an initial estimate vector 
to the nonlinear algebraic system before any iteration process is made. 
It is desirable that if the first estimate has failed to make the iteration 
convergent, another estimate is provided. In this subroutine, the 
initial estimate must be taken in the form of a second order polynomial. 
3. Function AUX 
This is a nonlinear algebraic equation system translated from 
equation (24), which has to be solved and which serves as an external 
input to subroutine ZSYSTM. 
14. Functions DIFF4 and DIFF8 
These are nonlinear algebraic equation systems translated 
from Eqns. 27 and 28, respectively, which have to be solved by the 
Finite Difference Method and which serve as external inputs to sub- 
routine DIFSOL. 
15. Subroutine LSQPL2 
This library subroutine, from the Naval Postgraduate School, 
Monterey, California, employs a least square best fit method to replace 
any function by a polynomial. 
ie oubroutine ZoYorM 
This library subroutine, from IMSL, solves nonlinear 
algebraic equations employing the Brown's method of iteration, which 


is a modified Newton's method and is developed in Reference 13. 
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Ge. LESTS BOR CONVERGENCE 

The test problems employed here are designed to verify the opera- 
tion of the computer program and test the convergence of the equation 
system. The two categories of test performed were linear and non- 
linear problems. 

An analytic solution and the Finite Difference Method solution were 
used to verify the Finite Element Method solution for linear and non- 
linear cases respectively. The following results (Tables 1, 2, 3) show 
that for an uniform beam of 100 inches in length, with a uniform load, 
the 8-element model gives sufficiently accurate results. These analyses 
show that as the number of elements is increased the finite element 
method results approach the theoretical results. The small difference 
between eight and sixteen element results (less than 1%) suggests the 
solution has converged sufficiently for engineering purposes. For 
symmetric problems (i.e., symmetric EI, k, s and q), advantage 
may be taken of symmetry. In this case, only half the system need be 
considered, thus greatly reducing computer effort. This reduction 
is especially beneficial for nonlinear problems. The time consumption 
may increase from 10 to 30 fold in going from a 4-~element model to an 
8-element model, depending on how good the initial a is. A brief 


description of the test program is shown in Figure 8. 
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Figure 8. Flow chart of a test program 
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Figure 10. Deflection of a beam ona sea floor model foundation 


compared to that of Winkler foundation. 


(Maximum difference is about 6%) 
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Table 1. Convergence of the Linear System 
(EIV") + KV =aq) 
Uniform Simply Supported Beam with Uniform Load 
Lb = 100 inte Er = 320 ibe ing, 
: q = 1 1b/in, k = 1 lb/in4, 
Deflection Given in Inches (BL = 0. 96) 


Location Finite Element Method Solution Theoretical 

x/L Number of Elements Number of Elements Solution 
8 16 

0.0 0.0 0.0 0.0 
0.125 0.016296 0.016310 0.016301 
0.250 0.029897 0.029919 0.029906 
On575 0.038839 0.038863 0.038880 
0.500 0.41950 0.041975 0.042178 


Table 2. Convergence of the Nonlinear System 
(BIV") + kVo =q 
Uniform Simply Supported Beam with Uniform Load 
L = 100 in, EI = 3 x 10° Ib-in“, 
q = 0.01 lb/in, k = 1 Ib/in@ 
Deflection Given in Inches ( BL = 0. 96) 


Finite Element Method Solution: 


Location Number of Elements Number of Elements 
x/L 4 8 
0.0 0.0 0.0 
0.250 0.1451D-03 0.0456D-03 
0.500 0.1483D-03 | 7 ¥ 0.1485D-03 
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Table 3. Solution of (EIV")" 1 eva =q 
F.E.M. Compared to F.D.M. Uniform 
Simply Supported Beam with Uniform Load 
ie= HOOsin, El = 6 9t23 1 0e beaee 
q = 2.25 lb/in, Number of Element = 8 
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Location k = 100 pie = "30 
oo) Vocal Dad Ws ip Ties 
ITMAX = 3 ITMAX = 3 
0.0 0.0 0.0 
®. 125 0. 18576rbD_03 OF 188397 D-03 
02250 0.340909D-03 0.245395 ).03 
0.375 0.442958D-03 0.448565D-03 
0.500 0.478469D-03 0.484450D-03 
Location k = 500 BL = 1.20 
x/L EF, E. 4. ae De \ © 
ITMAX = 8 ITMAX = 3 
0.0 O20 0.0 
Oa 0.185776D-03 0.188392D-03 
0250 0. 340900D-03 0. 345365 8)2038 
Oe BES 0.442946D-03 074485520405 
0.500 Oe ersva'Eys ID0)s' 0. 484436D-03 
Location k = 1000 Giae= 1242 
o/b pie Wie F.D.M. 
ITMAX = 4 ITMAX = 3 
0.0 0.0 Oan0) 
O25 0. 185769D-03 07 188385D-05 
0.250 0.340889D-03 0.345373D203 
023545 0.442931D-03 0.448536D-03 
0.500 0.478440D-03 0.484419D-03 
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x/L 


O70 

0. cs 
O22 50 
0.315 
Ur oUe 


Location 
x/L 


OO 

0. 125 
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Omens 
0.500 


k = 2000 
ITMAX = 4 


0) 

5 isis (fey (Wisi) 5) 
.340865D-03 
.442900D-03 
. 478407 D-03 


CO CO @ © © 


k = 3000 
ITMAX = 4 


0) 

. T8s74ob=03 
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gL = 1.69 
be IDE IME 
ITMAX = 3 


v0) 

else an 2203 
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Dm THE FINITE DIFFERENCE METHOD 

Since there is no analytic solution to the nonlinear case, the Finite 
Difference Method was used in this thesis to verify the results of the 
Finite Element Method. The development is restricted to uniform EI 
and k. The finite difference approximations are developed in Ref. 15. 
The fourth order derivative has the central difference computational 
molecule as follows: 


4 2 
(AV) = 1) fv(j-2) - 4V(j-1)46V(j)-4V(j+1)+V(j42)] + 0(1%) 


dx* ° 34 





where tis an element length 
Consider the differential equation of an uniform simply supported 
beam with uniform load ona uniform nonlinear foundation: 
tt 2 
HIV +kV =q 


ik! Z 
or EIV + EIV) =q 


k 
TEN? | 
ed k 
Let V = EIV and (Ee = A , then 


the above equation may be rewritten as follows: 


at | 


ie Rea 


-qe 0 
A computational molecule model applied at node j yields: 


1 V*(j-2) - 4V*(j-1) + 6V'(j) - 4V GFL) 4V (j42) 
x4 
+A V"*(5) - q(j) = 0 re 


Or. 
all 


pepezye4y (j-1) 6V (ij) ae Gel te) 
(26) 


Ae a, 
BIN Spe it) = oS Gl) = © 
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where j ranges from 1 to 2 for 4-element model and from 1 to 4 for 
8-element model. 


A. FOUR-ELEMENT MODEL 


V(-1) = -V(1) V(0) V(1) V(2) wel) V(0) 
e-— — = 
> 

x 


pati) 4 Veet A Xue gx 


I 
© 


(27) 
Beir Gu2Ne AuseAeieaag x 50 


B. EIGHT-ELEMENT MODEL 


= (1) V(0) V(1) V(2) v3) V(4) V(3) V2) V(1) V(0) 
e- -- —-o—__o_____o_—_o_____e—_—__»-—______9—____#___—__® 


Xx 


5V(1) - 4 V(2) +. V(3) +A X V%(1) - gx* = 0 


4 
-4V(1) + 6V(2) - 4V(3) + V(4) +A xty2(2) _ qX =0 


4.2 4 ee) 
Mal) . 4Viz) + ¢(5) - 4vqay A XV (3) - qx = 0 


IZ) OV (Sy OVA) ee (4) aX = 0 | 
Notice that the superscript * has been omitted for simplicity; the 
result V(j) must be multiplied by Sy to get the actual deflection. 
The above development is restricted to a uniform beam on a founda- 
tion with uniform k. In this case the finite difference method has 
distinct advantages over the finite element eee In the case of 


variable EI, k, and s the finite difference method is not easily developed 


and the advantage shifts to the finite element method. 
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VI. ILLUSTRATIVE PROBLEMS AND RESULTS 


This chapter presents the results of two investigations. 

The first investigation attempts to compare the behavior of three 
foundation models, the Winkler, the modified Winkler and the mixed 
mode foundations. The comparison is limited to the case of a simply 
supported uniform beam on foundations with uniform k and s properties. 
Table 4 gives values for the parameters considered inthis study. The 
value of gL varied from 1. 62 to 3.98 i.e., beams of medium length. 
The results obtained apply only to the particular simply supported 
systems considered. Other problems may yield other results. 

The second investigation seeks to establish that the Finite Element 
Program developed in this thesis can solve problems with variable EI, 
k, s and q. As a sample problem, the case of a free-free beam with 
variable stiffness, variable foundation modulus, and variable load was 
considered. The @L for this illustrative problem was 3.5. 

A. UNIFORM SIMPLY SUPPORTED BEAMS WITH UNIFORM LOADS 

ON SOIL FOUNDATIONS 

Figures 9 and 10 present the deflection ofa simply supported beam 
on soil foundations which are either iniand soil or sea floor sediment. 


The range of parameters considered is shown in Table 4 below. 
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a Table 4. Information Data 













Beam 
Properties 


Loading 
Conditions 





Foundation Characteristics 






q 
(lb/in) 







k s inland 
isin | lb soil and 
sea floor 
Sa Guerra: ene Te 


500 |1000| chosen 


as 


3000 |6000} examples 
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Appendix B contained the tabulated results of Figures 9 and 10 as 
well as data on the two foundations, inland soil and sea sediment. The 
results show that the difference between the Winkler hypothesis and 
other hypotheses is negligible for simply supported beams with uniform 
properties and loads. This may not be true for beams with other 


boundary conditions. 


B. A FREE-~-FREE NON-UNIFORM BEAM PROBLEM 
To show how the Finite Element Program can be used for a general 


case we consider the following problem: 








2 5000. 
4000.1lb/in 2 
Ib/in 
2000 
Ib/in ; 
k= ky + ks vi 
400 
lb/in? 





BG 





A symmetric problem was considered only to minimize the 


computer effort. The result is given by Table 5 and Figure ll. 
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Figure 11. Deflection of a free- a 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


The computer program presented in this thesis provides an accurate 
and reliable means for solving a variety of one-dimensional problems of 


beams on nonlinear foundations. The use of this program and efforts to 


increase its versatility are encouraged. 


A. REMARKS ON THE PROGRAM 

The following observations have been made during the development 
of the program and the investigation of various problems 

l. Curve Fitting 

It was assumed that the reaction of the foundation is proportional 

to vP , where V is the deflection and p any positive real number. 
Since the Finite Element Method via the Galerkin procedure is not easily 
applied directly for the case when p is not an interger, vP has been 


Zz 
replaced by a least square fit polynomial of the form b, + byV +b,V . 


] 3 


The reaction of the foundation must be zero when there is no deflection, 
therefore the greatest error may come from the non-zero residual 


constant b, (Figure 12). Suppose the least square fit curve r= by + 


3 : 
boV Tt ee intercepts the curve Ba vP at the point V=V. If the 
deflection of the foundation falls in the vicinity of V, the result is highly 
accurate. The result is less accurate, if the deflection is much larger 


or smaller than V. One must estimate how large or small the deflection 


will be, in order to select a "best fit'' curve. 
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Figure 12. Accuracy of curve fitting 


2. Initial Estimate of Deflection 
Iterative procedures are among the more popular schemes for 
solving systems of nonlinear algebraic equations. Such methods require 
an initial estimate of the solution. In this work the initial estimate was 
taken in the form of a second order polynomial in the subroutine 
VGUESS. Not all initial estimates will yield a solution; therefore it is 
desirable that if one initial estimate fails, it be replaced by another 
one. This is accomplished in subroutine VGUESS which changes the 
scale of the initial estimate function. : - 
3. Computational Considerations 
im any Computer eMortimene are two imajor concerns merc 


space required for a program and the time necessary for solution. 
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For the linear Finite Element problem the storage area required is 

proportional to m X M, where m is the number of coordinates and M 

is the bandwidth of the system. A change in mor M yields a correspond- 

ing linear change in the size of the computer storage. In addition the 

time for solution (for a Gauss elimination type scheme) is of the order 

of mx Mo for a symmetric system. These considerations show the 

disadvantages which result when mor M increase, and are well known. 
In the case of nonlinear problems the effects of increasing m 

and/or M on computational effort are even more pronounced. For 


example the space requirement associated with the nonlinear terms is 





aa x M. Hence an increase in m from ee to 7 increases the storage 
: ; aly) ere 5 eee ; ; 
in the ratio (4) . This is é times greater than in the linear case. 


Moreover, the computational effort is greatly affected since a large 
number of terms greatly increases the number of iterations. For 
example in the solution of an 8-element model (i.e., m= 18, M = 6) 
non linear problem, the storage for the array is about 8000 bits (single 
precision) and the computer time was about 60 seconds. When the same 
problem was modeled by a 4~element model (i.e., m= 10, M = 6), the 
storage was about 2400 bits (single precision) and the computer time 


was about 6 seconds; a reduction in storage in the ratio of 3 to landa 


reduction in computational time inthe ratio of 10 to l. 
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B. REMARKS ON THE WINKLER HYPOTHESIS 

Let a = kvP be the total settlement curve, i.e. the deflection 
versus reaction curve of a real foundation, which can be obtained by an 
éxperiment [Ref. 30]. Let the straight line r, = kV, where k is the 
foundation modulus, represent the behavior of the Winkler foundation 
(Figure 13). Let V correspond to the point at which the two curves 7 


and r intersect. Consider the area under the curve Be = kvP 
n 


between 0 and Vg. This area is equal to: 





Avs 
0 0 +] 
| rdv = kv? av = _* Vo" (30) 
ptl 
0 0 
Consider the functional: 
L Tf 1? 
+] 

F = | { (EIV) + (3 V) +—K_ vP e avytax (31) 

0 2 ptl 


where Lis the length of the beam. The general field equation (29) is 
the Euler's equation of the above functional. Therefore, to solve the 
equation (29) is the same as to minimize the functional (31). This 
functional represents the total potential energy of the system, including 
the beam, the foundation and the applied load. 

The third term under the integral is the energy due to the reaction 
of the foundation. According to Eqn. 30, the energy (31) is a function 
of the area under either the r, or r curve, depending on which founda- 


1 


tion is being used. Neglecting the effect due to the shear force in the 
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Figure 13. Comparison of a Winkler foundation 
to other types of foundations 
foundation for a while, for a specified beam and a given applied load, 
one can see that if the deflection is small, say less than V, the Winkler 
foundation deflects more than the nonlinear foundation does for equal 


supporting effort. Now, if the effect of the shear force in the foundation, 
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ee ee 


t 
ic. V a is taken into account, one can expect that the deflection 


= 
~O 
should be less. The shear force of the foundation depends on (sv) ¢ 


112 
which is usually much smaller than (EIV) and may be neglected. This 


discussion agrees with the results in two examples of two different 


types of soils given in Chapter VI, part A. 


Ce RECOMMENDED FUTURE STUDIES 
1. Following a similar procedure, the problem of plates or shells 
on two- or three- dimensional nonlinear foundations may be investigated. 
2. This thesis dealt with nonlinear foundations which are either 
continuous or discontinuous. Real foundations may be some combination 
of these foundations; hence there is a need to bridge the gap and consider 


foundations where the deformation is partly localized and partly 


continuous. 
3. Inthe present investigation the nonlinear term kv? was 
approximated by the second order polynomial, pet Dhl b,V". This 


restriction results from the finite element approximation of V(x) by 


m 
> V;G, (x). Future studies might consider ways to treat kVvP or in 
re 1 : : 

fact the more general case of an arbitrary nonlinear function of V, 


aarectly. 
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A. GENERAL FORMULATION 

This formulation is limited to plane stress theory only. Let u(x, y) 
and v (x, y) be the displacements at a certain point of the foundation, in 
the x- and y- directions.’ In general, an approximate solution can be 


“obtained by expanding u(x, y) and v(x, y) in a finite series: 


Ij 
u(x, y) >> U(x) oily) , 
i=1 


n (32) 
v(x, y) = » Nees) ~p Ay) ’ 
ce J J 
j=l 
where the dimensionless functions gy) and by) are known and 


the functions U, (x) and V j(x), which have the same dimensions as 
u(x, y) and v(x, y), are unknown. 
From the theory of elasticity, in the two-dimensional plane stress 


case, the stresses and strains are related as follows [Ref. 3]: 


Oxx : hen, a Veeyy_ ) | 
vee 
- Ve 
s E (33) 
2 : iniiien cia! Vecra) 
1 -% 
ae 
Txy = Tyx = Y xy 
2(itve) 
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é = oe 





xx ——<—<6— 
Ox 
wy ay 
: oy ee. 
Vxy = Yyx = oe + Ov 





Oy (x 


where E. and v¥; are the modulus of elasticity and the Poisson's ratio of 


the foundation, ae and ¢ are the strains in the x- and y- direction 


and Vey is the shearing strain. 


In terms of the series (32), the relations (33) may be rewritten: 








m nN, 
Er y a ! 
Thee iWeluee. @ eV. | 
2 : 1 j 
1. ? eg | _ J 
f j=l 
Oyy = Ef we. ! S ’ (35) 
-[ vp) v,9, 7, mj 
lw 421 j=! 


r mM n 

= T. = f ’ 

Poe = io oe 
2(1+ ve) Uzl jz! 

where the "prime" denotes the derivative of the function with respect 


to its own argument. 


From the foundation model shown in Figure Il, consider now an 


elementary strip of length Ax (Figure 14). The functions U, (x) and 


a) can be obtained from the equilibrium conditions. According to 
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— =< i. 


Figure 14. Force system ona strip 


of the foundation 
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1 the total work of all 


Lagrange's principle of virtual displacements, 
internal and external forces acting on this strip due to any virtual dis- 


placement is zero. Let: 


: T (x,y) = x Tey) , (36) 
{= 

v (x,y) = z me vy) (37) 
| Fe 


be the virtual displacements, where U, and Ve are arbitrary constants 
at that location x of the strip. 

Since there is no force applied on the surface of the foundation, 
and the bottom of the foundation rests ona rigid base, there is no work 
done on these two faces. It is assumed that all properties of the founda- 
tion remain constant in the z-direction, hence there is no pressure 
gradient in the z-direction and therefore the total work done on two 
faces perpendicular to the z-direction is zero. The system of forces 
which must be taken into account is shown in Figure 14. 

The external forces are the result of normal stresses Oxx? 


i + somes Ax, the shearing stresses yy, 


__ OTxy | Ax and the distributed forces, X (x,y) in the x-direction and 
Ox 


Y(x, y) in the Y-direction. 


The internal forces are the result: of normal stresses and 


oy 


shearing stresses Gey The work done by these internal forces is 


A virtual displacement is any displacement which is consistent 
with the forces and constraints [Ref. 2]. 
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the strain energy of the strip [Refs. 4, 9]. 


The total work done onthe strip by any m +n virtual displacements 


is then given by the following expression: 











= H H 
T —— 
| ( — Ax) (wdy) U, - [' ce 2 ) (Axwdy) 
0 cs 0 
H 
ui X (x,y) (Axwdy) U, + | etsy Ax) (wdy) v:; (38) 
0 0 Ox 
Hi | H 
= | ( Oyy Fics) Ax ewdy je, | Y (x,y) ( Ax wdy ) Me = 0 
0 0 
eet, 2,3,. » 1 ’ 7} Sehba2e-3s ,n) 
where 
<i = Uj io) age ee, 
yy a DS ise Ue 2 “(y) (39) 
oY oy Ox 
e. = ia - V. P (y) (40) 
yy 2 J 
y 
LZ 
(2 _ Ax) (wdy) u. = work of external forces ( OT xx Ax\(wdy); 
0 Ax Ox 
Hi ae 
j. ( Cxy vagy ( Ax wdy) = strain energy due to enea sens stress © CT xy) 


ie .. | 
pila Ax) (wdy) 4 = work of external forces (0 Txy Ax) (wdy) ; 
0 for x Ox 


( ‘©. ) (&x wdy) = strain energy due to normal stress 
{ i o 


dio) 


la 
x u. 
{ (x,yX Ax wdy) en | Y(x, y) ( Ax wdy) Vj = work of body forces. 
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After dividing through py wAx and substituting u., vV., V5 


j xy? €yy by 


(36), (37), (39) and (40), the expression (38) may be rewritten as follows: 


H ; H 


m lel 
: > ¢ es if 
>I | — Ui Yily) dy - ( exy 0, ~ ily) dy + JoX (x,y) Uyy)ay 
i=l 0 ox 0 


+ Y If Vx Ww dy - i OyyV. ply) d 
—— ') W ly )dy ; yyV; ¥ jly) dy) 








j-! 
i i (xy) iF Y j (y) dy}= (41) 
0 
or 
m =. Jel H H ( 
> We. i (0 Geer Pi (yondy - ‘ Txyy! (y) dy + \ X (x, y)P .(y)dy] 
; 0 “0x 0 0 J 
i 
n 7 | lel 
uy Doty, an DE xy y ly aye , Ce ay Soy ie 
i= “Ox 


(42) 


Since u and ME are arbitrary, the above expression leads to the two 


equations in the x- and y-directions: 


a Hi H 
{ "06 xx pily)dy - J xy @ (y)dy u i ey a ys Aye (43) 
0 x 0 0 


lel - 
i 0 Txy L -(y)dy - i Oyy Y ly )dy + S Y(x, y) Y (y)dy = 0 (44) 
—, ? ; j 0 j 
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Substitution of (35) into (41) and (42) leads toa system of mtn 


second order differential equations in U; (x) and V s(x) as follows: 





H 
{ X , dy = 045) 
0 


f M 
S 
§ 

+ 
i}: 
eae. 

~& 
— 
are 
Qu. 
<i 
+ 


O [ OL. @ ; + ve Y ; Ww, dy 
_ ) | ie ne: 
. Zinc) > as 


0 


H 
_ at [Vv Ube. DY, ] Y, dy + { YY, dy = 0 


These are the two governing equations of a two-dimensional 


foundation which has a depth H. 


we. oF ECIAL CASE 
For the limited scope of this thesis, let the horizontal displace- 
ments u be negligible. This means that any virtual displacements in 
the x-direction are also negligible. Since U are arbitrary constants, 
Eqn. 36 implies that 9%,(y) must be identically zero. Only the last 


equation of the system (45) remains and it reduces to: 
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m H 
H t 
2? f ( WON Baie eho see Vv : vay | 
% Ox 2(1+V;,) Jo ata ay... ; SG Yi Yj 





H 4 

- + Y (x, -dy =0 4 
j (py)  j y (46) 

where j= 1,2,3,...,n, and Ee and Ve are assumed to be constant. 


The above expression may be rewritten: 





H H 
Dit k( “s \ Yi; dy) V; lee 2s i 4 ws dy] ¥, 
2 | Ox 2(1+v¢) 0 l-v4 Yo wer 


lel 
+ \ Yer Y dy = 0 (47) 
J 
0 
or 
m 
! t 
r, : Vey V; ) E Kj Ne (4 8) 
1-1 
j= 1, Gaara 
where 
H 
0 
is the reaction of the foundation. 
E e | 
Sil = so { Vi Way (50) 
2(1+Ve¢) 0 : . 
H 
KS 5 = Ei ‘ ’ ‘ 
l- v¥ i wy ey (51) 
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The coefficients S. 


ie A lbs/in?, are the two characteristics 


lbs/in, and Ki j 


of a foundation [Ref. 5]. The former accounts for the shearing stress 
distribution in the foundation material, and the latter is equivalent to the 
Winkler constant of proportionality. They are however, dependent on 
how the p. functions are chosen and are consequently not independent 
of each other. 

For simplicity, assume that the foundation has an infinite depth, 
which is usually true for any soil layer compared to the deflection of 
most structural systems. Let the one-dimensional function w(y) 
satisfy the boundary conditions: 

w(o) = 1. (52) 
~(oo) = 0. (53) 
Then the displacement at any point of the foundation is: 

v(x, y) = V (x) wy) (ee) 
in which V(x) is the deflection at the surface of the foundation, where 
y = 0. Furthermore, assume that the displacements decrease ex- 
ponentially with the depth of the foundation, then Wy) may be 
expressed as: 


Soe 
Wiy)se ” (55) 


where A is a known constant which depends on the properties of the 


foundation. The values of s and k* are then given by (50) and (51): 


E o8) E 
se : ee = f lbs/;, (56) 
2(1+V,) 42(1+V;,) 
0 
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BE 1 
k* = f {- (ene ee = AE Ibs/, 3 (57) 
7 0 


2(1-v#) 


Equation (48) may be rewritten 
’ 
ee) ene (58) 


where s and k*are given by (56) and (57). 


Gr NUMERICAL APPLICATIONS 
In accordance with soil mechanics, taking the length of a loaded 
area as unity, Reference 7 gives an empirical relation between bearing 


load and settlement of a foundation as follows: 


V (59) 





where w is called the influence coefficient, which depends on the 
shape of loaded area and the properties of the foundation. The values 
of wp are givenin Table 6.[Ref. 6]. Comparison of the Winkler 
hypothesis 

P= kV (60) 
to equation (59) yields: 


(61) 


Eqns. 57 and 61 show the similarity of the theoretical analysis and 
the empirical result. The value A can be computed by equating 
fea to (61): 


h= 2H (62) 
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Table 6 
Influence Coefficients 


(After Z. Wilun and K. Starzewski [Ref. 7]) 


Flexible foundation Rigid foundation 


Shape of Settlement |Settlement [Mean 
of centre of corner Settlement a 
foundation 
Ho 
O. Pe point o 
Circular 1.00 {perimeter 


Square 


e 


Rectangle 


Ze | 
5x l 
10 x 1 


100 x 1 
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For the purpose of this thesis, consider a loaded area of shape 
wx Ax, where Axis taken as a unit length, one inch, and w is 
equal to the width of the foundation (Figure 14). Assume that w 
ranges between 5 inches and 10 inches. From Table 6, 1.83 < p< 2.25. 
mere! = 2c, then, from (62) XA = 4. In order to determine s and k, 
the modulus of elasticity and the Poisson's ratio of the foundation must 
be known. 

The modulus of elasticity of rocks ranges from 100 to 10,000 
eet mn/ im" = Lk1bf /in“) and the Poisson's ratio from 0.1 to 0.3 
(Refs. 10, 7]. Deep sea sediment cores have modulii of elasticity 
ranging from 0 to 12.4 psi[Ref. 11]. The modulii of elasticity of 
beach soils vary from 30 to 570 psi, while those of inland soils vary 
from 3,000 psi to 11,500 psi[Ref. 7]. Most soils have a Poisson's 
ratio near to 0.5 [Ref. 12]. The values of s and k computed from 
(56) and (57) for various types of soils are given in Table 7. 

Table 7 
s and k Values Depending on Various Types 


of Soils (Loaded Area of Shape on Ie 
Be ee ec) L) 


Sea Floor Beach Soils General 
Sediment Inland Soils 
Es, psi 12.4 S02570 | . 3000-11, 500 
Ve O25 OFS 45 
s, lbs 18.6 50 5, 000-21, 000 
k, psi os: 20S S60 1800-7200 
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APPENDIX B 


DATA AND RESULTS FOR PROBLEMS OF CHAPTER Vl ebeeeR Ts 


-_ 


The following example refers to Figure 23-10 of Reference 30. 
Assume that the piles used for this soil have 24-inch diameters. The 
reaction for the foundation is 

r= (Test load) w lbs/in 
Ome 
With D = 24 inches and w = 6.5 inches we have the following table. 
Table 8 
Deflection Versus Reaction ofa Natural Soil 


From Experimental Data 


[Ref. 30] lbs /in 
Deflection Test load 
(inches) (tons) 
OO ZO: 57406 
G25 40. 1149.2 
0.40 50; 1436.5 
0.60 60. lH24. 5 
he sl 70, ZO & 
120 1. . 2097. 6 


A curve of this data is shown in Figure 15. 
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r(kip/in) r(lbs /in) 





inland sea floor 
soil 7 model 
9. 
33 Px 
6. 
"(Sn 
VS 
3. 
A; Vi(in) 


C) Experimental data curve 


[] Least square fit 


() Assumed foundation behavior curve 
AN Approximate Winkler hypothesis, k = 3000. Deyn 


Figure 15. Foundation reaction versus deflection 


(After Spangler, M.G., [Ref. 30]) 
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A. A TYPICAL INLAND SOIL 

Figure 15 shows reactions versus displacements of a typical inland 
soil (read scale on left hand side). 

Assume that the Winkler foundation modulus is about 3000 eer 

The deflection results, obtained by the Finite Element program, of 
a simply supported uniform beam with uniform load resting on this type 
of soil are given by Table 9. 

Column 1: Solution of (EIV")_ +kV =q 

Column 2: Solution of (EIV")_ + r(V) = q where r (V) is replaced 
by a least square fit curve of the reaction versus deflection curve given 
foeoure 15, 


tT! 
Column 3: Solution of (EIV") - (sV ) + r(V) =q where s and k 


corresponding to a particular foundation is given by Table 7. 


B. AN ASSUMED SEA FLOOR MODEL 

Since there is insufficient experimental data of sea sediment, it 
will be assumed that the reaction versus deflection curve of sea 
sediment has the same shape as that of inland soil. According to 
Table 4 and Reference ll, assume that the Winkler foundation modulus 
of sea floor is about 8. 3 Tee. It will be assumed that the inland soil 
reaction versus deflection curve is scaled down by 8.3/3000 to obtain 
the sea floor reaction versus deflection curve. (Figure 15. read right 


hand side scale). 
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The deflection results of a simply supported uniform beam with 
uniform load resting on this type of sea foundation model is given by 


Table 10. 
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Table 9. Compared Foundation Deflection 
of Various Hypotheses for a Typical Inland Soil 
Represented in Figure 15 


fee a = O0lin,) El =) 30x 10° ipeease Ge= cece > lbs jin 


-_ 


ik 3 AHO Ths ine, 4 = OOO lbeie 


Location (2) Modified (3) Mixed mode 
x/L (1) Winkler Winkler foundation 
0.0 0.0 0.0 0.0 
0.125 0.003690 0.350190D-02 0. 349588D-02 
0.250 0.006390 0.603699D-02 0. 602663D-02 
Gms25 0.007898 0.744817D-02 0. 743554D-.02 
0.500 0.008370 0. 788845D-02 0.787515D-02 


2 
by L = 100 in, BI = 3 x10 lbs-in, q = 22.5 lbs/in 
2 
k = 3000 lbs/in , s = 6000 lbs/in 


Location (2) Modified (3) Mixed mode 
x/L (1) Winkler Winkler foundation 
0.0 0.0 0.0 0.0 
0.125 0.005513 0.524425D-02 0.522 12 D202 
0.250 0.007740 0.728582D-02 0.726411D-02 
0.375 0.007988 0.749146D-02 0. 746383)D-02 
0.500 0.008010 0.739400D-02 0. 739234D-02 


2 
sc) = 100 in, El=3.x 10" Ibs-in , q = 225 lbs/in 
k = 3000 Ibs/in”, s = 6000 lbs/in 


Location (2) Modified ; (3) Mixed mode 
x/L (1) Winkler Winkler foundation 
0.0 0.0 0.0 0.0 
0.125 0.055125 0.536074D-01 075336516 )-01 
0.250 0.077400 0.749189D-01 0.746865D-01 
Oe315 O07 92875 0.773826D-01 0.772876D-01 
0.500 0.080100 0.764925D-01 0. 764673D-01 
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Table 10. 


Z 
L =100in, FIlI=3x 10! lbs -in 


a) q = 1 lbs/in 


Location 
x/L 


0.0 

Oel25 
ez 5 0 
ms (5 
0.500 


Winkler 


0 

platol7 i220! 
.241054E-01 
Rol etan 220 | 
cv Ao 0c peo 


ooo oO O&O 


b) q = 2.25 lbs/in 


Location 
x/L 


Vea 

Qe2 5 
0.250 
Oreo 5 
0.500 


c)q= 4.5 


Location 
x/L 


O20 

eel Z5 
02.250 
Om375 


Winkler 


0 

Hovol36 6201 
2942371E-01 
#7403035E-0] 
5 PS olihew B10)! 


ooo OC O&O 


lbs/in 


Winkler 


O70 

Oa5 72276 -01 
0.108474 
0.140727 


Modified 
Winkler 


0 

. 130091D-01 
.238239D-01 
.309054D-_01 
. 333630D-01 


oo oO 2 © 


Modified 
Winkler 


. O 

VEG (ee sNbistv 
-538043D-01 
.698002D-01 
piooot bt 


o-oo "© © 


Modified 
Winkler 


0.0 

O25 9t542 5-01 
0.108346 

0. 140568 


v9 


Compared Foundation Deflection 
of Various Hypotheses for a Sea Floor Soil 
Represented by Figure 15 


Mixed mode 
foundation 


58 

, 1300299204 
.238124D-01 
- 308904D-01 
.333467D-01 


S250 Oo 


Mixed mode 
foundation 


0 

-293644D-01 
2517 oro 04 
.697660D-01 
~ 753150D-01 


ooo CO O&O 


Mixed mode 
foundation 


O70 
0.591253 b= Onl 
0.108292 
0.140498 





dbece=s0- 25 1bs/in 


Location 
ay Winkler 
0.0 0.0 
@125 0.8884616E-01 
0.250 ONG 27 12 
02375 0.211090 
0.500 O; 2eqece 


e) q = 9 lbs/in 


Location 
x/L Winkler 
O20 O20 
Onel25 0.118455 
0.250 0.216949 
Oe 15 0.281454 
0.500 0. 303843 


tf) q = 15.75 Ibs/in 


Location 
x/L Winkler 
0.0 0.0 
Orel2 5 0.207297 
07250 0.379660 
One 5 0.492545 
0.500 0.531724 


eq = 22.5 Ibs/in 


Location 

x/L Winkler 
0.0 0.0 
oe 25 0.296138 
0.250 0.542371 
e315 Oe vs005 
0.500 0.759606 


Modified 
Winkler 


me 
.893441D-01 
. 163656 
.212344 

. 229248 


o-Oo30 © © 


Modified 
Winkler 


0 

. 119966 
-219677 
.285174 
- 307884 


oO Oo O O 


Modified 
Winkler 


58, 

. 214630 
993 300 
i POET «7 
M55 1182 


oo Oo OC O 


Modified 
Winkler 


0 

Soha Sh 
215550 
. 747221 
. 806886 


oOo ©o @ eo 
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Mixed mode 
foundation 


a0) 
-892998D-01 
. 163574 
2122 3m 
SCALES), | Bie 


oO Go 2 OO @& 


Mixed mode 
foundation 


. 0 

. 119906 
en ora, 
. 285028 
307726 


=] ©: Oo © © 


Mixed mode 
foundation 


20 

oo, 
oo 9S Om 
. 510204 
. 550667 


Oo Oo @ OC © 


Mixed mode 
foundation 


. 0 

. 3 Poele 
2 (oe te 
. 746800 
. 806430 
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APPENDIX 


COMPUTER PROGRAM NOMENCLATURE 


A complete listing and description of all variables used in the 


program is not practical. The items listed in this appendix are common 


to several areas of the program and will assist the reader ina study of 


the program. For convenience items are listed by variable type. The 


definition or function and dimension, if applicable, of each item is 


given. 


Pe INTEGER CONSTANTS 


KA 


KEI 


KQ 


KS 


Pee oT 


KETEST 


KOQTEST 


ITMAX 


NBC 


NCODE 


maximum number of data of foundation modulus 
maximum number of data of beam flexural rigidity 
maximum number of loading conditions 

maximum number of various type of soils 

specified foundation modulus for test program only 
specified beam flexural rigidity for test program only 
specified loading condition for test program only 
maximum number of iterations : 
number of boundary conditions imposed 


code number, if NCODE is not equal to 1, the shear force 


in the foundation material is not taken into account 


8] 





NDOF number of degrees of freedom 


NELEM number of elements 
NELMAX maximum number of elements 
NNOD number of nodal points 

NSIG number of significant figures 


B. REAL CONSTANTS 


ALPHA foundation modulus k 

El beam flexural rigidity 

Q applied load 

SSOIL shear force ina soil foundation 
TL total length of the beam 

X beam element length 


ee ECTORS 

ASAV(10) a set of foundation modullii 

B(10) coefficients of a least square fit polynomial 
EF(10) element force vector 

EISAV(10) a set of beam flexural rigidities 

IDBC(20) identification number of bounda ry conditions 
PO(10) a set of power p's in the term kV | 
QSAV(10) a set of applied loads 


SK(10) a set of foundation modulii of natural soils 


55( 10) a set of shear forces of natural soils 
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memo (20) 
THETA(Z0) 
V(20) 


WA(600) 


D. MATRICES 
DC(4, 4) 

DM(4, 4) 

DN(4, 4, 4) 


EK(4, 4) 


ICORR(20, 4) 
SYSC(20, 20, 20) 


BS yoK(20, 20) 


system force vector 
slopes at the nodal points 
deflections at the nodal points 


working vector 


C3}. - matrix, referred to Eqns. 23 and 24 

M:, - matrix, referred to Eqns. 23 and 24 

Ni jk - matrix referred to Eqns. 23 and 24 
element stiffness matrix Key referred to Eqns. 23 
and 24 

correspondence table of local and global points 


system N; matrix, referred to Eqn. 24 


J 
system Ciz, Mj, and K;,_ combined matrix, 


referred to Eqn. 24. 
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